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control scheme. Moreover, we study the effect of the finite trion lifetime and of the imperfections 
of the unitary evolution such as off-resonant excitations and the nonadiabaticity of the driving. We 
calculate the total error of the operation on a spin-based qubit in an InAs/GaAs quantum dot system 
and discuss possible optimization against the different contributions. We indicate the parameters 
which allow for coherent control of the spin with a single qubit gate error as low as 10~*. 
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I. INTRODUCTION 

The spin degree of freedom of an excess electron in 
a charged semiconductor quantum dot (QD) has been 
proposed as an attractive candidate for the use as a 
qubiliii^. The advantage of storing the logical values in 
spin states is their long coherence time^. In addition, it is 
possible to optically induce charge dynamics dependent 
on the spin state of an electron via the Pauli exclusion 
principle and optical selection rules^i^ii. Coupling be- 
tween charge states via statioii^ and interband^i^ electric 
dipole moments allows one to perform quantum condi- 
tional operations on two spins. In this way, one can im- 
plement single- and two-qubit gates in QD systems. The 
utilization of the optical control methods leads to shorter 
switching times, on a picosecond time scale, in compari- 
son with nanosecond magnetic control of the spin, which 
is essential for the implementation of quantum informa- 
tion processing schemes. The capability of encoding and 
manipulating information at the single-spin level is of 
great importance and has been experimentally demon- 
strated recently: the generatioi*iS and optical controliiii^ 
of the spin coherence together with a possible read-out 
of the state of a single confined spin in a QD systemic 
can make the implementation of a quantum computer 
feasible. 

A promising scheme of quantum optical control of a 
spin in a single QD was recently proposed in Ref. [ij. It 
has been shown that coupling to a trion (charged exci- 
ton) state leads to an arbitrary rotation between the two 
Zeeman-split spin states. In this way, optical coherent 
control of a spin in a QD via adiabatic Raman transi- 
tions is possible. This control protocol does not require 
an auxiliary fourth state which was needed in a similar 
scheme previously proposedi^, removing the requirement 
for the transfer of the electron between two QDs and the 
delocalized hole state. 

In realistic experiments, the QDs used for the proposed 



implementation are embedded in a solid state matrix, and 
confined carriers interact with the phonon bath, leading 
to a loss of coherence. In optical spin control schemes, 
spin rotation is achieved by spin-dependent charge evo- 
lution. Therefore, the spin decoherence in these schemes 
results mainly from the lattice response to the evolving 
charge densityii^ii^ii^. This effect has been studied^^ for 
a control scheme using an auxiliary statei^ but not for 
the scheme of Ref. [ij that seems to be advantageous in 
some respects. Because of the different scheme of control 
the phonon-induced errors in the latter case have a dif- 
ferent form and can be attributed to two channels: pure 
dephasing and phonon-assisted trion generation. Fur- 
thermore, since in this control procedure the spin rota- 
tion involves a finite trion occupation (unlike in the other 
scheme^-), in addition to the phonon interaction the fi- 
nite trion lifetime can also lead to decoherence, because 
of the nonzero probability of the radiative decay of the 
trionl^. Moreover, the operation in the ideal case should 
be performed adiabatically and the imperfect adiabatic- 
ity of the evolution will contribute to the total error of 
the operation. Finally, since this control scheme involves 
spectral selection of transitions (in contrast to the other 
one), off-resonant terms can also lead to an unwanted 
leakage to the trion state and result in large discrepan- 
cies in the desired spin rotation. It is not clear in advance 
which of these factors (if any) will dominate the decoher- 
ence under specific driving conditions. So far, for the spin 
control scheme of Ref. Il4l . only the error resulting from 
imperfect adiabaticity of the driving and from the finite 
trion lifetime was evaluatediS, but the impact of phonons, 
off-resonant excitations, and the joint decoherence effect 
have not been studied. 

In this paper, we study the combined influence of the 
phonon and photon environments and of the imperfec- 
tions of the evolution on an optical spin control scheme 
based on an off-resonant coupling of the spin states to 
a trion state in a doped semiconductor QD^*^. As a 
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single-qubit gate, an optically induced arbitrary rotation 
of the electron spin state is considered. We show that 
the interactions with phonons and photons are the domi- 
nant sources of decoherence in this system. The phonon- 
induccd decoherence has two origins: pure dephasing and 
phonon- assisted trion generation, with a different depen- 
dence on the duration and the detuning of the optical 
control pulses. We show that by slowing down the evolu- 
tion one can considerably decrease only the error due to 
the pure dephasing, while the phonon-assisted trion gen- 
eration part still results in a large operation error. The 
contributions to the error resulting from the finite life- 
time of the trion and the imperfections of the evolution 
are also studied. We calculate the total error of the op- 
eration and study the nontrivial interplay of the different 
contributions. In particular, we show that for moderate 
detunings (~ 1 meV) and long pulse durations the er- 
ror is dominated by phonon-rclatcd effects. Wc indicate 
the optimal conditions for the qubit rotation for which 
the error is considerably small, with values even below 
10~*, which is essential for coherent quantum control, 
and discuss the possible optimization against particular 
contributions to the error. 

The paper is organized as follows. In Sec. HH we intro- 
duce the model for the qubit based on the spin states of 
a confined electron in a charged QD. Next, in Sec. IIIIl 
we describe a single-qubit rotation scheme. In Sec. IIVI 
we study the imperfections of the evolution. Section |V] 
describes the general perturbative method for describing 
the effects of the environment on an arbitrary operation 
on a logical qubit. Section IVll contains the results for the 
error contributions due to the interaction of the carriers 
with the phonon and photon environments. In Sec. IVIIl 
we calculate the total error of the spin rotation and dis- 
cuss the interplay and possible optimization of the dif- 
ferent sources of the error. Section IVIIII concludes the 
paper with final remarks. In addition, some technical 
details are presented in the Appendixes. 



II. MODEL SYSTEM 



where the first term (control Hamiltonian) describes the 
carriers and their interaction with the classical driving 
field (the laser beam) and the second is the sum of the 
free phonon iJph and photon i?rad contributions. The last 
part, V — i/c-ph + Hc-ind, describes the coupling of the 
carriers to the environment, where the first term denotes 
the carrier-phonon interaction and the second describes 
the coupling of the carriers to the photon field. 

The control Hamiltonian for this system in dipole and 
rotating wave approximations is 



He = ^e„|n)(n| 

71 

h 



^ n„(t)e'(""*+'^")(|0)-|-|l))(2|-|-H.c. 



.n=0,l 



where e„ are the energies of the corresponding states, 
LUn are the laser frequencies, and 7n are the phases of 
the pulses. The frequencies luq and uji have to satisfy 
the Raman conditions, namely that the detunings from 
the corresponding transition energies £2 — should be 
the same. To this end, we set ujn = (£2 ~ ^nj/f^ ~ A 
(n = 0, 1), where A is the common Raman detuning. We 
perform a unitary transformation to the rotating frame 
with |n) = e*(""*-'^''V) = 0>1)- We can set 70 = 7, 
7i = because only the relative phase is important. In 
the rotating frame, the Hamiltonian reads 



He = h^m\ + -n^{t) (e'^|0)(2| + e-^^|2)(6|) (1) 



-n,[t){m\ + \m)+Hi 



where 

H'c = Jr!i(Oe^^-*|6)(2| + ^f^o(t)e-^^-*+*^|i)(2| +H.C. 

contains oscillating (off-resonant) terms which can be 
treated as a perturbation to the ideal evolution gener- 
ated by the Hamiltonian given in Eq. ([T]). 



We consider a single semiconductor QD charged with 
one additional electron. A magnetic field applied in the 
X direction [see Fig. [^a)] causes a Zeeman splitting Ab 
between the electron states with spin-up and spin-down 
(with respect to the direction of the magnetic field) which 
define the two logical qubit states |0) and |1) in the pro- 
posed schem&ii. These two states are off-resonantly cou- 
pled (with a detuning A) to a trion state |2) by cr+- 
polarized laser pulses with real amplitudes ^o{t) and 
Vli(t) and with different phases [see Fig. [lja,b)]. These 
two electron spin states |0) and |1) together with the 
trion state |2) compose a three- level system, known in 
the literature as a A systerr^. 

The Hamiltonian of the system is given by 

H = He + i/onv + V, 




FIG. 1: (a) Schematic plot of the QD with the direction of 
incidence of the (T+-polarized laser beam and the orientation 
of the magnetic field, (b) Three-level system in a single doped 
QD. 
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The free phonon Hamiltonian has the form 

k 

where and f3k are phonon creation and annihilation 
operators, respectively, with corresponding frequencies 
tjfc, where k is the phonon wave number. The unper- 
turbed photon Hamiltonian in the absence of charges 
reads 



H, 



ad 



q,A 



with photon creation and annihilation operators cl^^ and 
CqA, respectively. Here, lu'^ — c|q|/ni. is the photon fre- 
quency, with the speed of light in vacuum c, the photon 
wave number q, and the refractive index of the semicon- 
ductor medium n^; X labels the polarization. The un- 
perturbed evolution of the system and of the decoupled 
environment is described by He + Hp^ + H^i^d ■ 

The interaction of the carriers with the environment in- 
cludes the coupling to phonons and photons. The carricr- 
phonon interaction reads 



ph 



E 



ki)(n'|^/„„'(fc)(/3fc+/9lfc) (2) 



with coupling constants /„„' (fc) having the symmetry 
fnn'{k) = /*/„(— fc). The states |0) and |1) correspond to 
a single electron confined in the same QD structure and 
differ only by the spin orientation. Therefore, they have 
the same orbital wave functions, and thus the coupling 
constants /oo(fc) and fii{k) are the same. These states 
have different spins so that foi{k) and /lo(fc) would de- 
scribe a 'direct' phonon-assisted spin flip, mediated by 
the spin-orbit coupling^^. However, in the optical spin- 
control schemes, the effect of this process is many orders 
of magnitude weaker than the decoherence due to the 
dynamical response to charge evolution^^. Therefore, we 
neglect this coupling and set these coefhcients equal to 
zero. 

Initially, before the arrival of the pulses, the lattice 
is in a configuration where one electron is present in 
the QD surrounded by a lattice deformation (a polaron- 
like stat o^^'^'^1^^ ). We redefine the phonon modes to ob- 
tain the ground state of the carrier-phonon system cor- 
responding to this new lattice equilibrium. In terms of 
the new phonon operators 



/oo(fc) 



the interaction with the lattice modes reads 
i7e-ph = \2){2\Y,F22{k){b^+bl^) 

k 

|i>(2|^Fi2(fc) (b^ + blA +H.C. 



(3) 



where F22(fc) = /22(fe) - /oo(fe) and Fi2{k) = 
/i2(fc)e~*("i*~'''f^ Additionally, there is a polaron-hke 
energy shift which is included in the control Hamilto- 
nian Hq. The intcrband off-diagonal phonon coupling 
[the second term in Eq. ([3])] has a negligible effect due 
to energetic reasonsAZ and can be disregarded. For spec- 
trally narrow pulses, which are needed for the adiabatic- 
ity of the procedure, only acoustic phonons contribute to 
the dephasing. For overlapping electron and hole wave 
functions, the excitation of a confined trion does not in- 
volve considerable charge redistribution and the effect of 
the piezoelectric coupling is very weak^^. Therefore, we 
consider only the interaction with longitudinal acoustic 
phonons via the deformation potential coupling. 

We assume for simplicity that the trion state is de- 
scribed by a product of electron and hole wave functions 
^c(h) {'"') which are the same as those corresponding to a 
single confined carrier. This is a reasonable approxima- 
tion in the strong confinement limit, where the Coulomb 
interaction is of minor influence on the wave functions 
and leads only to energy renormalization effects which 
are included in the transition energies. The coupling be- 
tween the trion and phonons [the first term in Eq. ^] 
has the form (see Appendix Ej) 



J^22(fc)-/22(fc)-/00(fe) 



/ hk 

2pVci 



{D,^Di,)T{k), (4) 



where p is the crystal density, V is the normalization 
volume of the phonon modes, c\ is the longitudinal speed 
of sound, and -Do(h) is the deformation potential constant 
for the electron (hole). The form factor J-{k) depends on 
the geometry of the wave functions ^'c(h) (*")• We assume 
Gaussian wave functions 



c(h)l 



exp 



(5) 



where Zo(h) is the confinement size for an electron (a hole) 
in the xy plane, Iz is the common confinement size along 
z, and r±_ and are the corresponding components of 
the position. Then, neglecting the small correction re- 
sulting from the difference between the electron and hole 
confinement sizes, one finds 

^(fc) = e-('=-'/2)^-(fc.'./2)^, (6) 

where P — {l^ + l'^/2 and k±_ and kz are the compo- 
nents of the wave vector in the xy plane and along z, 
respectively (see Appendix lAl for details). 

The carrier-photon interaction Hamiltonian in the ro- 
tating wave approximation reads (in the rotating frame) 



^-rad = ^^5qA4A[e'"°'|0)(2| 



(7) 



+e' 



q,A 

i(cJit-7) I 



|i)(2| +H.C. 
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Deformation potential coupling Do — -Dh 8 eV 



Crystal density p 5360 kg/m^ 

Speed of sound (longitudinal) ci 5150 m/s 

Wave function width in-plane I 5 nm 

Wave function width in z direction 1 nm 

Trion decay rate F 1 ns~^ 

Zeeman splitting /iAb 1 meV 



TABLE I: System parameters used in the calculations, 
with the coupling constants 

a—l ' 

Here, a denotes Cartesian components, da is the inter- 
band dipole moment, eg and Cj. = are the vacuum 
dielectric constant and the semiconductor relative dielec- 
tric constant, respectively, and e(^)(q) is the unit polar- 
ization vector. 

In Tab. HI the material parameters (corresponding to a 
self- assembled InAs/GaAs system) are given. 

III. UNPERTURBED SPIN ROTATION 

In this section, we present the formal description of 
the spin rotation procedure^ without the interaction 
with the environment. In the ideal case, the evolution 
is slow and may be described by invoking the adiabatic 
thcorcm^^ . 

In the considered three-level system [Fig. [Ijb)], it is 
possible to perform an arbitrary rotation of the electron 
spin via the intermediate trion state |2). To show this, 
one sets in the control Hamiltonian Hq [Eq. ([T])] 

A = e(t)cos[20(i)] 

and 

f7o(t) = f}(t)cos/3, ni{t) ^n{t)sinl3, 

where 

n{t) = e(t)sin[20(i)]. 

Here, 9(i) = ^^'^{t) + IS? is the (time-dependent) effec- 
tive Rabi frequency and 

sin^ </)(t) = i I 1 , ^ I . 

For A > 0, switching the pulses off corresponds to (j) ^ 0. 

To ensure an adiabatic evolution, the angle /3, defined 
via tan/3 = f2i/r2o, should vary slowly in time. We 
choose r^o s-nd l^i to have the same envelope shapes, so 
that /3 becomes time independent. 



We introduce new basis states 

\B) = e''^cos/3|6) +sin/3|i), 
\D) = e*'^sin/3|0) - cos/3|i), 

which are superpositions of the qubit states |0) and |1) se- 
lected by the laser pulses, where only the bright state \B) 
is coupled to the trion state and the orthogonal dark state 
\D) remains unaffected. The Hamiltonian He which gen- 
erates the ideal evolution [Eq. ([T])] now reads 

Hc = hA\2){2\ + ^nit)i\B){2\ + \2){B\) 

and has the instantaneous eigenstates 

\aoit)) = \D), (8a) 
|ai(i)) = cos</>(t)|B) -sin</>(i)|2), (8b) 
\a2{t)) = sin (l){t)\B) + cos (j){t)\2) (8c) 

with the corresponding eigenvalues 

\o{t) = 0, 

\i{t) = -he{t) sin^ (j){t) = J (a - v/A2 + f72(t)) , 

X2{t) = hQ{t) cos" c^{t) - ^ (a + ^/A^ + n^t)') . 

The system evolution is realized by the change of the so 
called tipping angle (j){t) [see Fig.[21[a)], i.e., by the change 
of the pulse amplitudes. The condition for adiabaticity is 
a slow change of (j){t) in comparison with the rate of the 
adiabatic motion given by the effective Rabi frequency, 
|0(i)| ^ If the adiabatic condition is met, the 

state of the system (initially a combination of |0) and |1)) 
remains in the subspace spanned by the two eigenstates 
|ao) and \ai) during the whole process. 

The evolution operator Uc{t) in the absence of the 
environment perturbation (in the basis \D), \B), |2)) has 
the form 



Uc{t)^\ e-'^i(*) cos (/)(*) e-*^=^(*)sin0(i) , 
\0 -e~*^iWsin(/)(t) 6"*^^ W cos 0(i) / 

(9) 

where 

1 /■* 

A«W = -r / drA„(T), (n= 1,2). 

From Eq. it is clear that after the operation, when 
the pulses are switched off {(j) 0), an arbitrary initial 
electron spin state will have acquired a phase in the \B) 
component with respect to the orthogonal dark superpo- 
sition \D). The resulting unitary transformation in the 
qubit subspace can be written as 

TT I \ -i-A^S-n -'^i T ■ • ^1 ^ -> 

Uc[oo) = e 2^ =cos— I-ism— cr-n. 
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FIG. 2: (a) The tipping angle and its derivative as a function 
of time for a detuning /lA = 1 meV and a pulse duration 
7-p = 5 ps. (b) The adjusted amplitude of the control pulse 
for the n/2 rotation about the z axis as a function of detuning 
for different fixed pulse durations. 



quirement of adiabatic evolution, as well as the effect of 
the oscillatory terms contained in H'q in Eq. ([T]) and ne- 
glected in the discussion presented in the previous section 
(the former has also been studied in Ref. [l9l ). 

A perfectly adiabatic evolution does not transfer the 
qubit states |0) and |1) into the trion state |2), which 
is only slightly occupied during the gating. In realistic 
experiments, the parameters cannot be changed infinitely 
slowly, so that there is a nonzero probability of a jump 
from the ideal instantaneous (adiabatic) state to one of 
the other states. Representing the exact system state in 
terms of the adiabatic eigenstates [Eqs. ([8al -c)]. 

|^)=^c„(i)e-'^"(*)|a„(0), 



where I is the unit operator and a is the vector of Pauli 
matrices in the original qubit basis |0), This transfor- 
mation corresponds to a rotation through an angle Ai(oo) 
about the axis n, 

n=[— cos7sin(2/3), sin7 sin(2/3), — cos(2/3)], 

which depends on the ratio of the pulse envelopes and on 
the relative phase of the pulses. 

For convenience, we write the general initial state of 
the spin qubit in the form 



cos -IB) 
2 ' ' 



sin -ID) 
2 ' ' 



(10) 



where d and ip are angles on a Bloch sphere. We will also 
need two states orthogonal to IV'o), which we choose in 
the form 

|V.i)=sin^|i?)-e»^cos^|i?), |^2) = |2). 

We assume that the rotation is performed using Gaus- 
sian control pulses 



17(i) =r!,exp( - — 



where $7* is the amplitude of the control pulse and Tp 
its duration. In the following discussion, we will treat 
the pulse duration Tp and the detuning A as tunable pa- 
rameters, while the pulse amplitude fi* will be adjusted 
to achieve the desired rotation of the qubit. The ampli- 
tude as a function of detuning for a 7r/2 rotation about 
the z axis [growth direction, see Fig. [T^a)] is plotted for 
different pulse durations in Fig.[2Ub). 



IV. IMPERFECTIONS OF THE UNITARY 
EVOLUTION 

Before we study the impact of the environment, let us 
discuss the limitations imposed on the driving by the re- 



one finds the equation for the probability amplitudes'^, 
Cra{t) = -^e'[^'"(*)-^"(*'l(a™(i)|a„(0>c„(0- 

n 

The state jag) is time independent and represents the 
dark state \D) decoupled from the trion state |2). More- 
over, (ao(i)|ai(i)) — 0. Therefore, to the leading order, 
the only unwanted transition is to the state |a2), which 
becomes the trion state after switching off the pulses. 
The corresponding amplitude is 



rfi e«[^=(*)-^i(*)l(a2(i)|ai(<))ci(0). 



If the initial state is |(^o(^)), as given by Eq. pU|) . then 
ci(0) = cos(i?/2) and 



rfi e^[A2W-Ai(t)]^(i), (11) 



The other source of imperfections in the controlled evo- 
lution are the rotating (off-resonant) terms contained in 
the Hamiltonian which reads in the new basis 



,i(ABt-7) 



sin (3 — e 



-i{A-Bt—y) 



cos^/3 



m2\ 



+ -n{t) cos(ABt - 7) sin(2/3)|B)(2| + H.c. 



We assume that the correction to the unitary evolution 
resulting from these terms is small and treat them pertur- 
batively. The effects of the additional Hamiltonian Hq 
may be of two kinds: additional unitary rotation within 
the computational space and leakage to the trion state. 
The former can be taken into account when designing 
the control pulses and compensated by a suitable mod- 
ification of the control parameters. Therefore, we treat 
only the latter as an error. The amplitude for the trion 
excitation is given by 
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„off 

C2 



dt{2\H^{t)\^/jo) 



1 ,q 

- cos - sin(2/3) / e*[^^(*)-^i(*)ll7(i) cos(ABi - 7) cos[20(i)] 



dte'^^^'^n{t) cos0(i) 



~i(ABt-7) 



„i(ABt-7) 
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FIG. 3: The error due to the unitary corrections as a function 
of detuning for two values of the Zeeman splitting Ab and 
three values of the pulse duration: Tp = 5 ps (solid lines), 
Tp — 10 ps (dashed lines), and Tp = 30 ps (dotted lines). 

where HQ{t) = UQ{t)HQU(.{t) is the relevant Hamilto- 
nian in the interaction picture with respect to the perfect 
evolution described by the evolution operator in Eq. 

The error due to both unitary corrections described 
above is given by the total probability of leakage to the 
trion state, 

X — I ^na I off 1 2 

and is plotted in Fig. [3] for a 7r/2 rotation about the z 
axis as a function of detuning for two values of the Zee- 
man splitting and different pulse durations. If the values 
of the detuning approach the Zeeman splitting Ab , the 
error becomes very large, since one spin state is then al- 
most resonantly coupled to the trion state. This induces 
a large trion occupation which inhibits a coherent spin 
rotation in this scheme. If the evolution is fast (short 
pulse durations Tp), the adiabatic condition is not met 
for small detunings, which results in larger errors. 

In addition, the unitary error shows oscillations, vis- 
ible in Figs. [3]Ja,b), which are due to the nonadiabatic 
contribution. In order to understand their origin, let us 
note that, in spite of the smooth Gaussian pulse enve- 
lope, the tipping angle evolves in a step- wise manner [see 
Fig. Ufa)] (especially for strong pulses). Thus, (j){t) has 
two peaks of opposite signs at t — ±ii, where ti is a cer- 
tain time, depending on the pulse duration Tp (ti « 10 ps 
in Fig. [2ja)]). Hence, according to Eq. pT|) . the proba- 
bility for a non-adiabatic jump is approximately propor- 
tional to sin^[A2(ti) - Ai(ti) - A2(-ii) -I- Ai(-ii)] and 
is therefore an oscillating function of the control param- 



eters. This oscillation of the transition probability re- 
flects the interference of the amplitudes for non-adiabatic 
jumps when the trion is switched on and off. 

For reasonably long pulses, one can indicate two 
regimes of the detuning where the operation on the qubit 
has a high fidelity. The values of the detuning have to be 
chosen either above the Zeeman splitting or below it and 
larger than 0.1 meV. To take advantage of the available 
fast optical control methods and insure the adiabatic- 
ity of the evolution, one can perform the operation with 
short pulse durations of several picoseconds and detun- 
ings of a few meV. 



V. ENVIRONMENT PERTURBATION DURING 
OPERATION ON THE QUBIT 

In this section, we summarize the general method for 
describing the effects of a perturbation due to the en- 
vironment on an arbitrary operation on a qubit. A full 
account of this approach is given in earlier worksiii^. 
Here, we give a brief overview for the sake of complete- 
ness and for reference in the following sections. 

The evolution in the absence of the perturbation due to 
the environment is described by the unitary evolution op- 
erator Uo{t) = C/c(i)'X'e"*^'=°'*. The effect of the interac- 
tion with the environment is calculated using the second- 
order Born expansion of the evolution equation for the 
density matrix. We include the effect of the driving field 
non-perturbatively and treat the carrier-environment in- 
teraction within a perturbation theory. 

The reduced density matrix of the qubit reads 

p{t)^Uo{t)[po + p^'HmUt), (12) 

where p^^) (^) jg ^j^e correction to the density matrix re- 
sulting from the interaction with the environment (in the 
interaction picture) and po the initial state of the qubit, 
which is assumed to be pure, po = |'i/'o)('0o|- The initial 
state of the system (qubit together with the environment) 
has the form po®pT, where pT is the thermal equilibrium 
state of the environment bath. 

To quantify the quality of the operation on a qubit, we 
use the fidelity^ 

F = (7/.o|C/t(oo)p((^)C/o(oo)|Vo)'/', (13) 
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which is a measure of the overlap between the ideal (pure) 
final state without perturbation, Uo{oo)\ipo) , and the ac- 
tual final state of the system given by the density matrix 
p(cxd). If the procedure is performed ideally, i.e. with- 
out discrepancies from the desired qubit operation, then 
F = 1. The fidelity loss 5 = 1 — is referred to as 
the error of the quantum gate. Inserting Eq. (|12p into 
Eq. ((T31), with po = |V'o)('0o|, one finds 



^ = -(V'o|p('noo)|Vo)- 



(14) 



The density matrix correction is calculated from a per- 
turbation expansion^^: 

P^'^i^)--^ fdr r dr'Tr^[V{T),[V{r'),po]], 

" Jto Jto 

where Tr^ is the trace with respect to the reservoir de- 
grees of freedom and V{t) = U^{t)VUoit) is the carrier- 
environment Hamiltonian in the interaction picture. It 
can always be written in the general form 



nn' : 



where Snn' acts on the carrier subsystem, Rnn' acts on 
the environment of phonons or photons, and Rnn' = 
Snn' = It is easy to see that Eqs. ([2]) and dH) 

have this structure. 

It is convenient to introduce two sets of spectral func- 
tions. The first is a family of spectral densities of the 
reservoir (phonons or photons), defined as 



Rnn', mm' {^) — J dt C^''^^ {Rnn' (t) Rmm') , (15) 

with the operator Rnn' transformed into the interaction 
picture Rnn'{t) = C/g (i)-R„„'?7o(0- The functions from 
the second set are nonlinear spectral characteristics of 
the driven evolution, 

Snn' , mm' \Y^,ni^Mm\Ymm'm, (16) 



where [ipi) span the subspace orthogonal to the initial 
state |?Ao) and 



Ynn'i^^) ^ dt Snn'{t)e 



(17) 



with Snn'{t) — U^{t)Snn'Uo{t). The various terms in 
Eq. describe transitions to different states orthogo- 
nal to the desired state {ipo)- In the long time limit for 
a time-independent system, all of them either vanish or 
turn into energy-conserving Dirac delta functions, restor- 
ing Fermi's golden rule for transition probabilities^^. In 
the general case, they are broadened due to time depen- 
dence. 



Using the definitions in Eqs. (|T5|) and ([T6)) . the error 
[Eq. (fT4| ] can be written in the form.— 



6 = 



nn' ,'mm' 



duj R 



nn' .mm' 



nn' ,mm 



(c). (18) 



The error can thus be expressed as an overlap between 
spectral functions, which are the "building blocks" for 
the calculation of the environment effects on the quantum 
evolution. A detailed derivation of Eq. (fT8|) can be found 
in Refs. [i3and[ll 

The perturbative approach described above obviously 
yields only an approximate description of decoherence. In 
the case of phonon-induced dephasing, comparisons with 
exact results (for ultrashort laser pulses)^ and with cor- 
relation expansion results^l show that the perturbative 
results are very accurate as long as the overall dephasing 
effect is weak: the inaccuracy is of the order of 6^ , where 
6 is some measure of decoherence, e.g., the fidelity loss, 
Eq. (fT8| . The same holds for the decoherence induced by 
the radiative decay of the trion, as we show in Appendix 

m 

Thus, to calculate the error of the quantum gate due to 
the interaction with the environment, one needs to derive 
the two spectral functions. These will yield a transparent 
spectral interpretation for the various contributions to 
the qubit dephasing and provide a possibility to seek the 
optimal conditions depending on the system properties. 



VI. DECOHERENCE MECHANISMS 

In this section, the different kinds of decoherence mech- 
anisms (the coupling of the carriers to phonons and pho- 
tons) arc studied. We apply the general theory intro- 
duced in the former section to calculate the total error 
of a spin rotation through an angle of 7r/2 about the z 
axis. The quantitative estimates are given for charged 
self- assembled InAs/GaAs quantum dot. 

A. Interaction with the phonons 

It follows from Eq. ([3|) that the carrier-phonon interac- 
tion is described by just one pair of operators, ^22 — |2)(2| 

and i?22 = Efc^22(fc)(&fc + ^Lfe). Therefore, one needs 
only two spectral functions: the spectral density of the 
phonon reservoir i?22,22(^) = ^ph(^) and the spectral 
characteristics of the driving 6*22, 22(1^) = '5'ph('^) to cal- 
culate the phonon-induced error 



Sph = J dLjRpi,{uj)Sph{uj). (19) 

Using Eq. (fT5|) . one finds the explicit form of the former: 
Rpbiio) = (20) 

^ [uBiiu) + l^22(fe)P [S{tO - LOk) + Situ + LOk)] , 
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FIG. 4: Spectral density of the phonon reservoir at diflterent 
temperatures T. 



where n-Q^uj) is the Bose distribution function. 

The function i?ph(w) (see Appendix [Al for a detailed 
derivation) is plotted in Fig. 2] for dilferent tempera- 
tures T. The negative frequency part of the function 
corresponds to phonon absorption by the carriers and is 
nonzero only for finite temperatures, while the positive 
part represents the emission processes and always has fi- 
nite values. The phonon spectral density has a cut-off at 
the frequency to ~ c\/l, which corresponds to the inverse 
of the time phonons need to traverse the quantum dot. 

The carrier part of the interaction Hamiltonian Sph = 
|2)(2| leads to a spectral characteristics containing two 
parts corresponding to the two orthogonal states: 



ph/ 



(21) 



= - sin^ -d 
4 



dt ( 



' sin^ (/)(t) 



^^g»c.tg-»[A2(t)-Ai(t)] gij^ [20(t)] 



1 , d 
- cos - 



The initial state, Eq. lfTU|). can be an arbitrary superposi- 
tion of the states |0) and |1), and the error depends on the 
choice of the initial qubit state. In order to obtain rep- 
resentative error evaluations, we average the error (thus, 
the contributing spectral function) over the angles (t?, Lp) 
on the Bloch sphere of the initial states, according to 

Sj(av)('^) = -r diSsini) dip Si(w). 

The averaged contributions to the spectral characteris- 
tics, corresponding to the two terms in Eq. (j2ip . read 



s^]^ ■.(lu) 

l(av) ^ 



ph 
®2(av) 



(^) 



dt e'"' sin^ (l){t) 

^^^^u.t^-^[^2{t}-A^{t)] J20(t)] 



Both contributions to the phonon-induccd error are in- 
dependent of p. 

We can derive approximate analytical formulas for the 
spectral characteristics under the condition that the con- 
trol pulses are much smaller than the detuning (il ^ A), 
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FIG. 5: (a),(b) Spectral function sj'^((x;). (c),(d) Spectral 
function (uj). In (b), (d) the line styles denote pulse dura- 
tions as defined in (d). 



which is met for detunings of several meV [see Fig. [l{b)] . 
Thus, in this detuning regime, we obtain 



(C) 
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exp ^ 
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|exp 
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2^/3A- 



■ exp 



The symmetric function Si^^{uj) [see Fig.[5l^a,b)] is cen- 
tered at a; = and, for longer pulse durations, covers the 
low frequency part (broadening l/^p). The broaden- 
ing of this function is independent of the detuning, and, 
for a fixed Tp, its area decreases for larger detunings. 
This part of the spectral characteristics corresponds to 
pure dephasing effects^^'^'^i^'^. The resulting error 
[Eq. will decrease for longer pulse durations as well 
as for larger detunings. 

The second part of the spectral characteristics {oj) 
[see Fig.[51[c,d)] is centered at u; w — A, and its area grows 
with time. Thus, the corresponding error contribution is 
proportional to the spectral density of the phonon reser- 
voir around the frequency corresponding to the detuning 
of the laser frequency from the trion transition, with some 
broadening due to time dependence. Moreover, the error 
increases for longer operations. Therefore, this contribu- 
tion may be interpreted as a real transition and describes 
the error due to phonon-assisted trion generation. 

In order to see this more directly, let us note that the 
function S2^{oj) is relatively strongly localized around 
LU = —A, compared to the range of variation of R(uj). 
Therefore, the corresponding integral in Eq. (|19p may be 
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approximated by its Markovian limit, 

The area of the spectral function 52'^'^) appearing in 
this formula may be calculated by noting that for any 
function h(t), one has 



duj 



dt e"^*/i(i) = 

duj J dtj dt'e"^^'-'"^h*{t)h{t') 
= 27r / dt\h{t)\^ . (23) 



Thus, we find 

= 2tt j dt i?pii(-A)i cos^ ^ sin^ [2(j){t)] . 

The expression under the integral is exactly the Fermi's 
golden rule formula for the probability that a transition 
from the state |ai) to \a2) will take place during the time 
dt due to the (diagonal) carrier-phonon coupling given 
in Eq. ([3])- Since the state [02) becomes the trion state 
after the laser pulse is switched off, this process indeed 
represents a phonon-assisted transition to the trion state. 

The resulting phonon-induced errors, averaged over all 
initial states, as functions of the detuning and pulse du- 
ration are shown in Fig. [51 The error due to pure de- 
phasing [corresponding to s5'^(cj)] favors longer pulse 
durations and larger detunings and strongly depends on 
the temperature [Fig. [6Ua,b)]. To perform the operation 
with an error smaller than 10~^, which allows for coher- 
ent quantum operation on a qubit, one needs a detuning 
larger than 0.06 meV (at T = 5 K). It is possible to avoid 
the pure dephasing error even for a fast evolution realized 
by pulse durations of several picoseconds and a detuning 
of 1 meV. 

The contribution to the total error related to the 
phonon-assisted trion generation 5^ [resulting from 
S2*^(w)] has a different behavior in comparison with the 
previous one [Fig. [ni[c,d)] and depends even stronger on 
the temperature. At low temperatures [T < 1 K), it 
decreases with growing detuning, but for higher tem- 
peratures (T = 1, 5 and 10 K), it initially grows with 
detuning and later decreases for considerably large de- 
tunings. The maximum values correspond to the situa- 
tion when the spectral function 82^^ (uj) covers the maxi- 
mum of the phonon density, and the error vanishes if it 
lies beyond the phonon density cut-off. The pulse du- 
ration dependence of this error differs from the one for 
the pure dephasing error. For relatively large detunings 
(hA > 1 meV), this error favors shorter pulse durations, 
which is typical for real transitions. 

In order to properly choose the conditions for the spin 
rotation that lead to a high fidelity of the operation, one 
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FIG. 6: (a) The dependence of the error due to pure de- 
phasing for growing detuning for Tp = 10 ps. (b) The pure 
dephasing error as a function of pulse duration at T = 5 K 
and for five values of detuning: hA — 0.066 meV (solid line), 
hA = 0.66 meV (longer dashed line), hA — 1.32 meV (shorter 
dashed line), hA = 2 meV (dotted line), and hA = 2.64 meV 
(dashed-dotted line), (c), (d) As in (a), (b) but for the error 
due to phonon-assisted trion generation, (e), (f) As in (a), 
(b) but for the total error resulting from the interaction of 
the carriers with the phonon reservoir. In (a), (c), and (e), 
the line styles denote temperatures as defined in (a). In (b), 
(d), and (f), the lines styles denote detunings as defined in 
(b). 



needs to take into account these two sources of error re- 
sulting from the carrier coupling to the phonon reservoir. 
The sum of these two errors is plotted in Fig. [6l(e,f). To 
achieve values of the error below 10""* for Tp — 10 ps, 
detunings of several meV are needed. For a small detun- 
ing (hA — 0.066 meV), the error decreases with growing 
pulse duration, but for larger detunings {hA = 0.66 meV 
and hA = 1.32 meV), shorter pulse durations are more 
favorable. Choosing the detunings larger than several 
meV can suppress the influence of the phonon environ- 
ment (hA > 4 meV at T = 5 K). 

It is also possible to excite the system above the tran- 
sition energies by choosing a negative detuning. In this 
case, the phonon-induced total error is larger in com- 
parison with the case for positive detuning (see Fig. [7]). 
The corresponding spectral characteristics (uj) is now 
centered at the positive frequency part of the phonon 
density, which has larger values than the negative part, 
finite even at T = K (see Fig. HJ. This corresponds 
to emission of a phonon, which is strongly temperature 
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FIG. 7: The phonon-induced error for negative detunings (a) 
at different temperatures for growing detuning and (b) for 
growing pulse duration. 



dependent and possible at zero temperature. Especially 
at low temperatures, the phonon-induced error is up to 
four orders of magnitude larger than the one for positive 
detunings. To suppress the influence of the phonon reser- 
voir, one needs larger detunings. The dependence of this 
error on the pulse duration also differs. Only for small de- 
tunings {hA < 1 meV), it is advantageous to use longer 
pulses. For relatively large detunings {HA > 1 meV), 
short pulses are favorable. 



B. Interaction with the photon field 

Since the considered procedure for spin rotation re- 
quires a trion occupation during the evolution, the radia- 
tive decay of the trion will result in an additional error. 
We calculate the effect of the photon interaction in the 
same manner as for the phonons including perturbatively 
the carrier-photon interaction. 

The relevant photon energies correspond to the semi- 
conductor band gap which is very large compared to the 
thermal energy. Therefore, one can use the zero tempera- 
ture approximation. The carrier-photon interaction con- 
tains the following operators acting on the carrier subsys- 
tem: ^02 = |0)(2|e^-°*, ^20 = ^0^2, ^12 = |l)(2|e^("i*-'^), 
and 5*21 = 5'j2- AH the resulting contributions, calcu- 
lated according to the general procedure in Sec. |Vl can 
be combined into a single spectral characteristics 



5.adH = 5]|(V'o|i;i,(c.)iv^.) 

i 



where 

yrad(t^) 



1 

V2 



dte 



—iujt Tj] 



ul{t) 



X [(cos/3|B) -Ksin/3|L>))e*"'°* 

-f (sin/3|B) -cos/3|Z?))e*"i*] {2\Uo{t) 



and the sum is taken over all states orthogonal to the 
initial state |V'o) [Eq. pU)) ]. 

The spectral density of the photon reservoir reads^S 



^-d(c^) = ' ' 3 , > 0, 



(24) 



and the contribution to the error due to the photon in- 
teraction has the form 



5rad = y 0^tJi?rad(w)S'rad(^^)• 



Using the definition ([T7]) and the explicit form of the 
evolution operator ©J one finds 



{MYL{U^M^ 



where 



- - sin i9 suv[2(t){t)] (e-"^°* cos (3 -I- e"*"i* sin (3) 



„j[v+Ai(t)] ^ 



i9 



cos^ -sin0(i)(e'"^i*cos/3- e-*""'sin/?) 



and 



s\^''{t) = -e-^[^^«-^iWlcos|sin2 0(i) 
x(e~*'^''*cos/3-He-"^i*sin/3). 



The spectral function S'i.ad(^) is centered at the laser 
frequency lo = wo,i ~ 1 eV, and its width is of the 
order of 1 meV or less. The photon spectral density 
-Rrad(i^) [Eq. (EH)] is very broad and may be assumed 
constant in the area of the overlap with S-[ad{^)- There- 
fore, we use the Markovian approximation and obtain 

^rad = i?rad(wo) / da;S'rad(w) with i?rad(t^o) = r/(27r), 

where F is the trion decay rate « 1 ns~^. The frequency 
integral can be performed using Eq. (j23p . Upon averag- 
ing over the initial states, one obtains the resulting error 
due to the carrier-photon interaction: 
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-^rad 



r / dt 



^ sm^[2(j>{t)] + i sin^ {1 + cos [(^i - uJo)t] sin(2/3)} 

+ - sin^ (/)(t) {1 - cos [{uji ~ ujo)t] sin(2/3)} 
3 



(25) 



10' 



10"- 



10 



...(a) T =5ps — 
- ■-■>... ^ 10 ps — 

30 ps ---- 


(b) M = 0.0066 meV ... 






;r--'---"*^0.666 meV '. 


"^0.66 meV 
6.6 meVs^ 



0.01 0.1 1 

hA (meV) 



10 

^r, (ps) 



100 



FIG. 8: The dependence of the radiative error for growing (a) 
detuning and (b) pulse duration. 



For large enough detunings (fi <C A), we can again 
obtain an approximate equation for the resulting error: 



-^rad 



V 8A2 



8\/2A4 



exp 



exp 



4 P 



1 



p(^i -^^o)^ 



The resulting radiative error [Eq. (|25p ] as a function 
of detuning and pulse duration is plotted in Fig. [8](a,b). 
This error decreases with growing detuning since the 
trion occupation is reduced. For small detunings, the 
trion occupation is relatively large since the system is ex- 
cited near the resonance, and the resulting error is grow- 
ing with pulse duration due to the growing probability of 
the radiative decay of the trion. In this regime, the error 
is linear in the pulse duration, (5rad ~ TpT. The contribu- 
tion to the error due to the finite trion lifetime depends 
strongly on the pulse duration only for relatively small 
detunings. For detunings of several meV, this error is 
constant in time. To insure a small radiative error, one 
has to properly choose a large detuning, while the pulse 
duration can be arbitrarily short. 



VII. INTERPLAY OF THE DIFFERENT KINDS 
OF ERRORS 

In the previous sections, we studied the detuning and 
pulse duration dependence of the contributions to the 
total error of the spin-based quantum gate. In this sec- 
tion, we calculate the resulting total error and discuss 
the interplay of and possible optimization against the 
constituent sources of the error. 
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FIG. 9: Interplay of the error contributions as a function of 

(a) , (c) detuning and (b), (d) pulse duration at T = 5 K. 
Different lines represent various error sources, as defined in 

(b) . 



We start the discussion with the dependence of the 
contributions to the total error of the considered spin ro- 
tation on the detuning for a fixed temperature (T = 5 K) 
and two pulse durations (rp = 30 ps and Tp = 5 ps) 
[see Fig. |9l[a,c)]. For small detunings, the trion occu- 
pation is large, and the dominant source of the error is 
the radiative decay of the trion. This contribution de- 
creases with growing detuning, and the phonon-induced 
error becomes the most important source of dephasing. 
If the values of the detuning approach the Zeeman split- 
ting, one spin state is almost resonantly coupled to the 
trion state, and the probability of the leakage to the trion 
state is high, which leads to a large error due to the uni- 
tary corrections, especially in the case of a short pulse 
(rp = 5 ps). In the detuning regime between 0.09 meV 
and 2 meV, the error resulting from the phonon coupling 
becomes dominant and is one order of magnitude larger 
than the radiative error. For detunings of several meV 
above the value of the cut-off of the phonon error, the 
only contribution which inhibits the coherent spin rota- 
tion is the trion radiative decay and has values between 
10-3 and 10-4. To achieve errors smaller than 10 , 
large detuning of several or tens of meV are needed. In 
this detuning regime, the limitation results from the fact 
that one cannot choose an arbitrarily large detuning, in 
particular not with a frequency corresponding to that of 
optical phonons, assumed to be well off-resonant in this 
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paper. Thus, the interplay of the different kinds of error 
leads to a nontrivial detuning dependence of the total 
error, which in general is dominated by the error due to 
the interaction with the phonon and photon reservoirs 

Let us now discuss the dependence of the particular 
error contributions on the pulse duration for a fixed tem- 
perature (T = 5 K) [see Fig. [Hl^bjd)]. For a small detun- 
ing {hA — 0.33 meV) and fast driving fields (short pulse 
durations of several picoseconds), the adiabatic condi- 
tion is not fulfilled, and the probability of the leakage 
to the trion state becomes very high, which results in 
large errors. To minimize the influence of the errors due 
to the unitary corrections, one has to choose pulse du- 
rations of at least a few picoseconds. In this regime, 
the total error is limited by the phonon-induced contri- 
bution. The second dominant source of the error is the 
radiative decay of the trion. Moreover, these two contri- 
butions are almost constant in time. For a large detun- 
ing {hA — 6.6 meV), the evolution is perfectly adiabatic, 
and the error due to unitary corrections does not appear. 
The phonon-induced error is very small and vanishes for 
pulse durations of a few picoseconds, since the detuning 
is far above the value of the cut-off of the phonon spectral 
density. The only important contribution to the total er- 
ror is due to the finite trion lifetime which is constant 
in pulse duration. One can see that it is advantageous 
to excite with larger detunings, but the variation of the 
pulse duration is not of great importance as soon as the 
adiabatic condition is met and the detuning is not close 
to the Zeeman splitting. 

In order to summarize the study of the various sources 
of decoherence, we calculated the total error of the spin- 
based qubit rotation for growing detuning [Fig. [TUTa')] 
and pulse duration [Fig. fTUTbl]. For detunings smaller 
than 1.3 meV, the error strongly depends on the pulse du- 
ration and is relatively large. Thus, to perform a rotation 
with a high fidelity, detunings of several meV are needed. 
Furthermore, for such detunings, the error is constant 
for growing pulse duration, which opens the possibility 
to perform the qubit rotation with pulse durations three 
orders of magnitude smaller than the lifetime of the trion 
(r^^ w 1 ns). For hA ~ 5 meV, the error may be as low 
as 10^** and is independent of the pulse duration. 



VIII. CONCLUSIONS 

We have studied a theoretical proposal of the opti- 
cal control of a single spin in a single doped quantum 
dot^'*. We have investigated the sources of error of a 
quantum operation on a spin-based qubit and have given 
quantitative estimates for the implementation of spin ro- 
tation through an angle of 7r/2 about the z axis in a 
self-assembled InAs/GaAs system. The dephasing mech- 
anisms resulting from the interaction of the carriers with 
phonon and photon reservoirs as well as the imperfections 
of the unitary evolution have been considered. 

We have shown that the interplay of the constituent 
sources of the error leads to a nontrivial dependence of 
the total error on the detuning, which is in general dom- 
inated by the errors due to the coupling of the carriers 
to the phonon and photon environments. Furthermore, 
small detunings or detunings approaching the Zeeman 
splitting should be avoided since they lead to large trion 
occupations which suppress an adiabatic and coherent 
control. Taking into account all contributions to the to- 
tal error, we showed that errors as low as 10^^ can be 
achieved for large detunings (~ 5 meV), while the pulse 
durations can in principle be arbitrary (but at least of a 
few picoseconds). 

Finally, it should be noted that the calculations were 
performed using simple Gaussian pulses with the pulse 
duration and intensity as the only tunable parameters. 
Further reduction of the errors is very likely to be possible 
with pulse optimization32i2^i^. 
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APPENDIX A: PHONON COUPLINGS, 
SPECTRAL DENSITY 

In this Appendix, we derive the effective coupling el- 
ement for the carrier-phonon interaction F22{k) as well 
as the resulting spectral density of the phonon reservoir 
i?ph(w) for the studied QD system. 

The general interaction Hamiltonian for confined 
states (restricted to the ground states of the carriers) 
reads 

^^int = y^y^ \DcJ^c{k)al^„aa,a - Dh^h{k)al^ah,a 

k a 

where cie(h),CT) ^l(h) <y ^^'^ annihilation and creation op- 
erators, respectively, for an electron or a hole in the con- 
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fined ground state with spin cr, and the form factors are 
given by 

/ + 00 
-OO 

(Al) 

For Gaussian wave functions as in Eq. ^ , one exphcitly 
finds by a simple integration 

JPo(h)(fc) = exp (^-^fci/2^h) - I'^lll^ ■ 

For a single electron state |0) = a^||vac) or |1) — 
al J vac) (I vac) is the empty dot state), one immediately 



finds 



(0|i/i„t|0) = (l|i/i„t|l) 



/ hk 

2pVc{ 



^e(fe) ibk + bl^ 



For the trion state \2) — a^, t^a\^ ||vac), one has 
(2|ffi„t|2) = 

E J ^ [2^e^c(fc) - DMk)] {b^ + bl, 



In order to further simplify the calculations, we ne- 
glect the difference between the localization sizes of the 
electron and hole wave functions. Taking into account 
the small variation of the electron and hole confinement 
widths Iq and Ih leads only to inessential quantitative 
corrections^. This leads to the form factors Tc{k) — 
J^ h(fc) = T{k) [Eq. ©] and to Eq . © with /oo = /ii = 

^D,T{k) and h2 - /^(2i?c - D^)T{k), and 
thus to Eq. (g]). 

With the isotropic acoustic phonon dispersion, we ob- 
tain the spectral density of the phonon reservoir R-p]^{Lo) 

[Eq. my- 

i?ph(w) = [uBiuj) + l] (-2^)3;^2 '^^ j ^^dCcosC 

X f dkk^ \F22ik)f [SiiU - LJk) + S{iJ + UJk)] , 



where the angles ij and ( denote the orientation of the k 
vector. We can rewrite it in the form 



where 



Rq 



h{D, - D^Y 
87r2pcf 



and the function g{Lo) is defined as 



-it/2 



d( cos ( exp 



-tt/2 



f 2. ■ 2. 

-^(^cos^C+^sm^C 



APPENDIX B: LINDBLAD MASTER EQUATION 
FOR THE TRION RECOMBINATION CHANNEL 



In this Appendix, we derive the results for the trion 
recombination channel in the Lindblad formalism and 
compare them with those calculated by means of the per- 
turbative method discussed in this paper. 

From the interaction Hamiltonian -ffc-iad [Eq. (7)], we 
derive the Lindblad equation^^ in the form 



p = T[ <T-pa+ - - ^pa+a- ] - i [i/ad,p] 



where 



-iujQt 



{\2){B\ + \2){D\)+e-^-^'{\2){B\-\2){D\)] 



and i?ad = itfc{t)UQ{t) is the Hamiltonian generating 
the adiabatic evolution, with Uc given by Eq. This 
equation is consistent with the perturbative approxima- 
tion in the sense that the latter is reproduced upon trans- 
forming to the interaction picture and performing an ex- 
pansion in carrier-phonon coupling. 



10" 



10" 



10"^ 
10"^ 

0.01 0.1 1 

hA (meV) 

FIG. 11: The radiative error calculated by means of the Lind- 
blad (solid lines) and perturbative (dashed lines) methods for 
growing detuning. 

The results from the Lindblad equation together with 
those calculated with the perturative theory are plot- 
ted in Fig. 11 for the initial state ji/'o) = |^) for the 
7r/2 rotation about the z axis. For small detunings 
[h/S. < 0.1 meV), where the error is relatively large, the 
perturbative method yields slightly larger errors, while 
in the case of larger detunings the results are the same. 
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